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ABSTRACT 

In  this  paper  we  study  a  history-boundary  value  problem  for  a  nonlinear 
conservation  with  fading  memory  in  one  space  dimension .  The  motivation  for 
studying  this  problem  is  an  earlier  work  by  C.  M.  Dafermos  and  the  author 
concerning  the  motion  of  a  nonlinear,  one- dimensional  viscoelastic  body. 

Using  a  variant  of  an  energy  method  applied  to  the  viscoelastic  problem  it  is 
shown  that  under  physically  reasonable  assumptions  the  nonlinear  conservation 
law  has  a  unique,  classical  solution  (global  in  time),  provided  the  data  are 
sufficiently  smooth  and  "small"  in  a  suitable  norm;  moreover,  the  solution  and 
its  first  order  derivatives  decay  to  zero  as  t  +  ".  The  proof  illustrates 
the  versatility  of  the  energy  method  combined  with  frequency  domain  techniques 
for  Volterra  operators. 

A  preliminary  analysis  based  on  current  work  of  R.  Malek-Madani  and  the 
author  is  presented  concerning  the  development  of  singularities  in  smooth 
solutions  of  the  conservation  law  (in  finite  time)  for  sufficiently  "large" 
smooth  data;  under  special  assumptions  it  is  shown  that  such  singularities 
necessarily  develop.  The  hope  is  to  apply  such  a  procedure  to  the 
viscoelastic  problem. 
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SIGNIFICANCE  AND  EXPLANATION 


Problems  arising  in  continuum  mechanics  can  often  be  modeled  by  quasilinear 
hyperbolic  systems  in  which  the  characteristic  speeds  are  not  necessarily  constant* 
Such  systems  have  the  property  that  waves  may  be  amplified  and  solutions  that  were 
initially  smooth  may  develop  discontinuities  ("shocks")  in  finite  time.  Of 
particular  interest  are  situations  in  which  the  destabilizing  mechanism  arising  from 
nonlinear  effects  can  coexist  and  compete  with  dissipative  effects. 

In  certain  cases  dissipation  is  so  powerful  that  waves  cannot  break  and 
solutions  remain  smooth  for  all  time .  A  more  interesting  situation  arises  when  the 
amplification  and  decay  mechanisms  are  nearly  balanced  so  that  the  outcome  of  their 
confrontation  cannot  be  predicted  at  the  outset .  A  dimensional  analysis  indicates 
that  the  breaking  of  waves  develops  on  a  time  scale  inversely  proportional  to  wave 
amplitude,  whereas  dissipation  proceeds  at  a  rate  roughly  independent  of  amplitude. 
Therefore,  when  the  initial  data  are  small  it  might  be  expected  that  the  dissipation 
effects  would  prevail  and  waves  would  not  break.  Results  of  this  type  have  been 
obtained  by  Nishida  [17]  for  a  model  problem  concerning  a  quasilinear  second  order 
wave  equation  in  one  space  dimension.  Nishida's  analysis  uses  the  method  of  Riemann 
invariants  and  is  therefore  restricted  to  one  space  dimension.  Using  energy  methods, 
Matsumura  [16]  has  studied  the  case  of  more  than  one  space  dimension,  and  was  able  to 
prove  the  existence  of  smooth  solutions  for  all  time  for  quasilinear  hyperbolic 
systems  with  frictional  damping.  The  necessity  of  the  presence  of  some  form  of 
frictional  damping  to  avoid  the  formation  of  "shocks"  in  finite  time  follows  from  a 
fundamental  result  of  Lax  [10]  .  An  assumption  of  this  theory  is  that  the 
constitutive  relation  characterizing  the  nonlinear  partial  differential  operator  is 
convex,  but  in  the  case  of  non-convex  functions,  as  may  arise  in  problems  in 
nonlinear  elasticity,  similar  'blow-up*  results  have  been  obtained  by  MacCamy  and 
Mizel  [12] j  general  results  for  nonlinear  hyperbolic  equations  have  been  obtained  by 
John  [8]  and  Kleinerman  and  Majda  [9] . 

A  different  and  subtler  dissipative  mechanism  is  induced  by  memory  effects  of 
elastico-viscous  materials.  Dafermos  and  Nohel  [5]  have  recently  developed  and  ana¬ 
lyzed  a  one-dimensional  nonlinear  model  for  the  homogeneous  extension  of  an  elastico- 
viscous  rod  whose  ends  are  free  of  traction.  This  equation,  which  is  a  quasilinear 
hyperbolic  integrodif ferential  equation  of  Volterra  type,  is  based  on  the  assumption 
that  the  stress  is  a  particular  nonlinear  functional  of  the  strain  involving  two 
assigned  constitutive  relations,  and  an  assigned  stress-strain  relaxation  function. 
Shock  solutions  are  known  to  occur  for  this  model  in  the  special  case  in  which  the 
stress-strain  relaxation  function  is  identically  constant.  But,  by  combining  energy 
methods  with  frequency-domain  techniques  for  nonlinear  Volterra  equations,  it  has 
been  shown  in  [5]  that  the  nonlinear  integro-dif ferential  equation  which  describes 
the  motion  has  a  unique,  smooth  solution  for  all  time  provided  the  given  data 
(history  of  motion  and  external  body  force)  are  sufficiently  smooth  and  "small."  The 
decay  properties  of  the  solution  and  of  its  derivatives  have  also  been  analyzed.  The 
hypotheses  for  this  theory  include  all  currently  used  constitutive  and  stress-strain 
relaxation  functions.  This  work  generalizes  studies  by  MacCamy  [11]  who  used  the 
method  of  Riemann  invariants,  and  by  Dafermos  and  Nohel  [4]  and  Staffans  [20]  who 
used  the  energy  method?  in  these  the  stress-strain  functional  involved  a  single 
assigned  constitutive  function.  The  energy  method  developed  for  these  problems  was 
used  recently  by  Slemrod  [18]  to  study  the  evolution  of  smooth  solutions  of  a 
mathematical  model  for  thermoelasticity  where  the  dissipation  arose  through  the 
influence  of  thermal  damping. 

A  natural  and  difficult  question  for  the  viscoelastic  problems  is:  do  shock 
solutions  develop  (i.e.  do  the  waves  break  in  finite  time)  as  the  smooth  data  become 
sufficiently  large?  This  question  was  answered  in  the  affirmative  by  Slemrod  [19] 


J 


for  the  aforementioned  thermoelastic  problem.  For  the  more  subtle  viscoelastic  model 
the  problem  is  far  from  solved  in  its  full  generality. 


The  purpose  of  this  paper  is  to  study  a  simpler  model  problem  consisting  of  a 
nonlinear  conservation  law  with  memory  ((1.1)  below),  together  with  a  prescribed 
history  function  for  t  <  0,  and  a  prescribed  boundary  condition.  This  problem, 
which  is  of  independent  interest,  is  of  first  order  (while  the  viscoelastic  problem 
is  of  second  order);  however  it  displays  the  crucial  features  of  the  more  complex 
problem.  Namely,  it  is  guasilinear  hyperbolic  and  in  the  absence  of  the  memory  term 
the  model  problem  reduces  to  Burger's  equation  which  motivated  the  above  mentioned 
theory  of  Lax  on  the  formation  of  shocks  in  finite  time . 

In  this  paper  we  modify  the  energy  method  developed  by  Dafermos  and  the  author 
[5] ,  and  obtain  the  global  existence  of  smooth  solution  for  sufficiently  smooth  and 
"small"  data  of  the  model  conservation  law  with  memory. 

He  then  formulate  an  approach  to  studying  the  formation  of  singularities 
(shocks)  in  smooth  solutions  in  finite  time;  an  advantage  of  this  approach  is  that 
whenever  it  is  applicable,  physically  meaningful  entropy  conditions  are  satisfied  by 
the  shock  solutions .  He  present  one  such  result  covering  a  wide  class  of  problems . 
The  ultimate  objective  is  to  extend  these  methods  to  systems  of  conservation  laws 
with  memory,  and  eventually  to  the  viscoelastic  problem. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC ,  and  not  with  the  author  of  this  report. 


A  NONLINEAR  CONSERVATION  LAW  WITH  MEMORY 


J  i  A  •  Nohel 

1.  Introduction .  In  this  paper  we  study  the  model  nonlinear  Volterra  functional 
differential  equation  (with  infinite  memory) 

t 

(1.1)  u*.  +  <p(u)  +  /  a'  (t  -  t)<i<u(t,x))  dr  -  f(t,x)  <-*  <  t  <  •,  0  <  x  <  l )  , 

c  x  -  X 

where  :  R  ♦  R  are  given  smooth  constitutive  functions,  a  :  to,")  +  R  is  a  given 

memory  kernel,  and  f  :  R  *  [0,11  ♦  R  is  a  given  function  representing  an  external  forcei 
subscripts  denote  partial  derivatives  and  '  _  d/dt .  The  motivation  end  the  assumptions 
under  which  (1.1)  is  studied  are  provided  by  the  more  complex  physical  problem  of  the 
extension  of  a  finite,  homogeneous,  elastoviscous  body  moving  under  the  action  an  assigned 
body  force.  The  viscoelastic  problem,  formulated  in  Section  2,  was  recently  studied  by 
Dafermos  and  Nohel  (51 j  references  to  earlier  literature  are  given  in  Section  2. 

The  model  problem  (1.1),  which  is  of  independent  Interest,  is  simpler  in  that  it  is  of 
first  order,  while  the  equation  of  motion  (2.8)  below,  is  of  seoond  order;  otherwise  (1.1) 
incorporates  the  Interesting  features  of  (2.8).  The  most  important  of  these  is  the 
following.  If  ♦  =  0,  t  ;  0,  (1.1)  reduces  to  Burger's  equation  (conservation  law  of  gas 

dynamics) : 

(1 .2)  ut  +  lp*u)x  “  0  ' 

note  that  (1.2)  is  quasilinear  and  hyperbolic.  It  is  classical,  see  Lax  [10],  that  the 
Cauchy  problem  consisting  of  Burger's  equation  and  the  initial  condition  u(0,x)  •  uQ(x), 
x  e  R,  does  not,  in  general,  possess  a  classical  smooth  solution,  no  matter  how  smooth 
(and  "small")  one  takes  the  Initial  datum  u 0;  if  <fi  is  convex  I*x  [10]  has  shown  that 
the  solution  develops  a  singularity  (shock)  in  finite  time  due  to  the  crossing  of 
characteristics.  More  precisely,  if  uj(x)  >  0,  uQ  smooth,  and  l  is  convex,  (1.2)  has  a 
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the  characteristics 


smooth  solution  for  all  t  >  0,  while  if  u^(Xq)  <  0  for  some  xQ, 
of  (1.2)  will  cross  and  shocks  will  develop  in  finite  time;  a  similar  result  also  holds  for 
systems  of  conservation  laws  in  one  space  dimension.  If  <f  is  not  convex,  e.g.  in 
nonlinear  elasticity,  similar  results  have  been  established  by  MacCamy  and  Misel  [12);  for 
recent  closely  related  literature  see  also  general  results  by  John  [8] ,  Kleinerman  and 
Majda  [9],  Nalek-Hadani  (13). 

The  purpose  of  this  paper  is>  First,  in  Section  2  we  formulate  the  problem  of  motion 
of  a  nonlinear  viscoelastic  body  as  analysed  recently  by  Dafermos  and  Nohel  [5]  using  a 
combination  of  energy  methods  and  properties  of  strongly  positive  kernels .  second,  under 
assumptions  motivated  by  the  viscoelastic  problem,  mm  show  in  Sections  3  and  4  that 
equation  (1.1),  together  with  an  assigned  periodic  boundary  oondition  and  an  assigned 
smooth  history 

u(t,x)  -  v(t  ,x) ,  -“  <  t  <  0,  0  <  x  <  1  , 

which  satisfies  (1.1)  for  t  <  0,  possesses  a  unique  classical  solution  for  all  t  >  0, 
provided  the  history  function  v  and  the  forcing  term  f  are  sufficiently  smooth  and 
"small*  in  suitable  norms.  This  result,  in  which  the  same  strategy  as  In  [5]  is  used, 
exhibits  the  dissipative  character  of  the  Integral  term  In  (1.1);  its  proof  serves  to 
illustrate  a  general  energy  technique  for  hyperbolic,  nonlinear  Volterra  problems. 

Finally,  in  Section  5  we  formulate  the  problem  of  development  of  singularities  In  finite 
time  and  mm  present  a  recent  result,  analogous  to  the  result  for  Burger's  equation,  in  an 
important  special  case,  (15) ;  the  proof  will  appear  elsewhere.  The  fact  that  the  question 
of  development  of  singularities  for  the  physically  interesting  viscoelastic  problem  remains 
open,  provides  the  principal  motivation  for  studying  the  simpler  model  problem  (1.1).  It 
should  be  noted  that  the  problem  discussed  in  Section  5  is  different  from  the  study  of  weak 
solutions  for  the  Riemann  problem  for  (1.1)  (Greenberg,  Ling  Hsiao,  and  MacCamy  [61, 


Dafermos  and  Ling  Hsiao  [3] ) 


2.  A  Hyperbolic,  Nonlinear  Volterra  Equation  In  Viscoelasticity .  Problems  arising  in 
continuum  mechanics  can  often  be  modeled  by  quasilinear  hyperbolic  systems  in  which  the 
characteristic  speeds  are  not  necessarily  constant.  Such  systems  have  the  property  that 
waves  may  be  amplified  and  solutions  that  were  initially  smooth  may  develop  discontinuities 
(“shocks")  in  finite  time.  Of  particular  interest  are  situations  in  which  the 
destabilizing  mechanism  arising  from  nonlinear  effects  can  coexist  and  oompete  with 
dissipative  effects. 

In  certain  cases  dissipation  is  so  powerful  that  waves  cannot  break  and  solutions 
remain  smooth  for  all  time.  A  more  interesting  situation  arises  when  the  amplification  and 
decay  mechanisms  are  nearly  balanced  so  that  the  outcome  of  their  confrontation  cannot  be 
predicted  at  the  outset.  A  dimensional  analysis  indicates  that  the  breaking  of  waves 
develops  on  a  time  scale  inversely  proportional  to  wave  amplitude,  whereas  dissipation 
proceeds  at  a  rate  roughly  independent  of  amplitude.  Therefore,  when  the  initial  data  are 
small  it  might  be  expected  that  the  dissipation  effects  would  prevail  and  waves  would  not 
break.  Results  of  this  type  have  been  obtained  by  Nishida  [17]  for  a  model  problem 
concerning  a  quasilinear  second  order  wave  equation  in  one  space  dimension.  Nishida' s 
analysis  uses  the  method  of  Riemann  invariants  and  is  therefore  restricted  to  one  space 
dimension.  Using  energy  methods,  Matsumura  (16)  has  studied  the  case  of  more  than  one 
space  dimension,  and  was  able  to  prove  the  existence  of  smooth  solutions  for  all  time  for 
quasilinear  hyperbolic  systems  with  frictional  damping.  Burger's  equation  (1.2)  shows  the 
necessity  of  the  presence  of  some  form  of  frictional  damping  to  avoid  the  formation  of 
"shocks".  The  rather  delicate  situation  of  thermal  damping  in  thermoelasticity  is 
discussed  by  Slemrod  [18]  . 

A  different  and  subtler  dissipative  mechanism  is  induced  by  memory  effects  of 
elastico-viscous  materials.  Dafermos  and  Nohel  (5)  have  recently  developed  and  analyzed  a 
one-dimensional  nonlinear  model  for  the  homogeneous  extension  of  an  elastico-viscous  rod 
whose  ends  a.e  free  of  traction.  Their  simple,  one- dimensional,  model  corresponds  to  the 
following  constitutive  relation,  suggested  by  the  theory  developed  ty  Coleman  and  Gurtin 
M! , 
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(2.1) 


t 

<J(t,x)  -  *(e(t,x>)  +  /  a'(t  -  T)iKe(T,x))dT  , 

where  cr  is  the  stress,  e  the  strain,  a  the  relaxation  function  with  '  -  d/dt,  and 
•P.  ♦  assigned  constitutive  functions.  The  relaxation  function  is  normalized  so  that 
aC*)  *  0.  When  the  reference  configuration  is  a  natural  state,  l(0)  »  $(0)  »  0. 
Experience  indicates  that  Vie),  We),  as  well  as  the  equilibrium  stress 

(2.2)  X(e)  d2f  We)  -  a(0)*(e) 

are  increasing  functions  of  e,  at  least  near  equilibrium  (|e|  small).  Moreover,  the 
effect  of  viscosity  is  dissipative.  To  express  mathematically  the  above  physical 
requirements,  we  impose  upon  a(t),  l(e),  We)  and  x(e)  the  following  assumptions: 

(2.3)  a(t)  e  W^'^tO,**),  a(t)  is  strongly  positive  definite  on  (0,*)  i 

(2.4)  We)  6  C3(-«,«),  WO)  -  0,  •  (0)  >  0  > 

(2.5)  ♦(e)  6  C2(-»,«),  W0)  -  0,  ♦■(0)  >  0  ; 

(2.6)  X’<0)  -  ^>*  (0)  -  a(0)W(0)  >  0  . 

Assumption  (2.3),  which  requires  that  alt)  -  a  exp(-t)  be  a  positive  definite  kernel  on 
[0 ,”)  for  some  u  >  0,  expresses  the  dissipative  character  of  viscosity.  Smooth, 
integrable,  nonincreasing,  convex  relaxation  functions,  e.g. , 

K 

(2.7)  a(t)  -  l  v^expt-^t),  >  0,  Wk  >  0  , 

k-1 

which  are  commonly  employed  in  the  applications  of  the  theory  of  viscoelasticity,  satisfy 
(2,3). 

We  now  consider  a  homogeneous,  one-dimensional  body  (string  or  bar)  with  reference 
configuration  (0,1]  of  density  P  -  1  (for  simplicity)  and  constitutive  relation  (2.1) 
which  is  moving  under  the  action  of  an  assigned  body  force  g(t,x>,  -<*  <  t  <  ",  0  <  x  <  1, 
with  the  aids  of  the  rod  free  of  traction.  We  let  u(t,x)  denote  the  displacement  of 
particle  x  at  time  t  in  which  case  the  strain  is  e(t,x)  »  ux(t,x).  Thus  the  equation 
of  motion  Pu^  a  +  Pg  here  takes  the  form  of  the  nonlinear  (hyperbolic)  Volterra 
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functional  differential  equation 


t 

(2.8)  Uf.f.  =■  V(u  >  +  /  a*  (t  -  T)i|i(u  )  dT  +  g,  -»<t<»,  0<x<1. 

X  _  XX 

The  physical  problem  of  the  motion  of  a  viscoelastic  body  suggests  that  the  history  of  the 
motion  of  the  body  up  to  time  t  «  0  is  assumed  known,  i.e., 

(2.9)  u(t,x)  -  v(t,x),  “®  <  t  <  0,  0  <  x  <  1  , 

where  v(t,x)  is  a  given  sufficiently  smooth  function  which  satisfies  equation  (2.8)  for 
t  <  0,  together  with  appropriate  boundary  conditions.  In  order  bo  show  that  the  motion  of 
the  viscoelastic  bar  remains  smooth  for  all  t  >  0,  the  mathematical  task  is  to  determine 
a  smooth  extension  u(t,x)  of  v(t,x)  on  (-“,*)  x  [0,1)  which  satisfies  (2.8)  together 
with  assigned  boundary  conditions,  for  -*  <  t  <  *. 

Dafermos  and  Nohel  [5,  Theorem  1.1]  establish  such  a  global  result  for  the  problem 
(2.8),  (2.9)  together  with  homogeneous  Neumann  boundary  conditions 

(2.10)  ux(t,0)  «  ux(t,1)  -  0,  -»  <  t  <  «  » 

these  are  shown  to  be  equivalent  to  the  statement:  the  boundary  of  the  body  is  free  of 
traction  (o(t,0)  ■  o(t,1)  -  0,  <  t  <  ",  0  given  by  (1.1)).  Their  global  result  is 

valid  for  sufficiently  smooth  and  "small"  external  body  forces  g  and  history  functions 
v.  Other  boundary  conditions  and  various  generalizations  are  also  considered. 

The  general  strategy  used  in  [5]  is  as  follows.  First  one  establishes  the  existence 
of  a  unique  smooth  local  solution  u  defined  on  a  maximal  Interval  (-®,TQ)  *  [0,1],  with 
the  property  that  when  Tg  <  m  a  certain  norm  of  the  solution  becomes  infinite  as  t  *  Tg; 
this  is  done  by  a  fixed  point  argument  (combined  with  a  standard  energy  method  for  linear 
problems)  on  a  suitably  chosen  abstract  space  of  functions.  Second,  energy  methods  are 
combined  with  properties  of  strongly  positive  kernels  to  show  that  due  to  the  viscous 
dissipation  of  the  integral  term  in  (2.8),  the  aforementioned  norm  remains  uniformly 
bounded  on  the  maximal  Interval,  provided  the  data  g  and  v  are  sufficiently  smooth  and 
small.  By  standard  theory  for  nonlinear  problems  this  means  that  Tg  “  +®  and  the  smooth 
solution  exists  globally  in  t.  This  part  of  the  analysis  involves  obtaining  a  priori 
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estimates  of  certain  norms  of  the  derivatives  (in  one  space  dimension,  up  to  and  including 
order  3)  directly  from  the  equation  (2.8),  It  is  here  that  it  becomes  convenient  to  use 
the  equivalent  form 

t 

(2.11)  u*,.  -  X(u  )  +  /  a(t  -  T)i|i(u  )  dT  +  g,  -*•  <  t  <  “,  0  <  x  <  1  , 

XX  X  Tx 

of  equation  (2.8)  (equation  (2.11)  is  obtained  from  (2.8)  by  integrating  by  parts  with 
respect  to  r  and  by  using  the  definition  of  the  equilibrium  stress  X);  assumption  (2.6) 
plays  a  crucial  role. 

For  the  special  case  i|>(e)  =  i(e)  various  global  existence  results  for  (2.8),  were 
established  by  MacCamy  [11],  Dafermos  and  Nohel  [4]  and  Staffans  [20].  The  assumption 
’ll  =  V  allows  one  to  invert  the  linear  Volterra  integral  operator  on  the  right-hand  side 
of  (2.8)  and  thus  express  ^  ( ux )  x  in  terms  of  ut^  -  g  through  an  inverse  Volterra 
integral  operator  using  the  resolvent  kernel  associated  with  a'.  One  may  then  transfer 
time  derivatives  from  ufct  to  the  resolvent  kernel  via  integration  by  parts.  This 
procedure,  which  was  also  discussed  by  Dafermos  [2]  in  a  recent  expository  paper  in  a 
number  of  special  cases,  reveals  the  instantaneous  character  of  dissipation  and,  at  the 
same  time,  renders  the  memory  term  linear  and  milder,  thus  simplifying  the  analysis 
considerably.  On  the  other  hand,  the  above  approach  is  somewhat  artificial:  By  inverting 
the  right-hand  side  of  (2.8),  one  loses  sight  of  the  original  equation  and  of  the  physical 
interpretation  of  the  derived  a  priori  estimates.  More  importantly,  the  physical 
appropriateness  of  the  restriction  V  m  is  by  no  means  clear. 

Remark i  The  present  normalization  of  the  kernel  a  with  a(")  “0  i3  different  from  that 
in  the  existing  literature  (see  [4],  [11],  [20]).  The  reader  should  note  a',  not  a, 
enters  the  constitutive  relation  (2.1)  as  well  as  the  equation  of  motion  (2.8).  In  the 
earlier  literature  in  which  only  the  special  case  ’ll  =  was  studied,  the  normalization 

a(t)  -  8.  *  A(t>  0  <  t  <  »  , 

2  1 

a(0)  «  1,  a—  >  0,  A  e  W  '  (0,~) ,  A  strongly  positive  was  used.  The  normalization  used 
here  is  crucial  for  generating  the  a  priori  estimates  directly  from  equation  (2.  Ill 
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(equivalent  to  (2.8)).  The  reader  should  note  that  the  present  normalization  and  (2.6) 
imply  that  0  <  a(0)  <  1,  if 

The  question  of  the  development  of  singularities  of  solutions  of  (2.8)  in  finite  time 
for  sufficiently  "large",  smooth  data  (which  have  been  observed  for  viscoelastic  bodies 
[1])  is  under  active  study.  Some  partial  results  with  $  5  ^  in  (2.8)  have  been  obtained 
by  Hattori  (7)  ?  for  a  viscoelastic  fluid  (see  Slemrod  [19]).  However,  in  the  general  case 
of  (2.8)  the  problem  is  far  from  settled.  For  this  reason  we  believe  that  the  approach  in 
Section  5  for  the  considerably  simpler  problem  (1.1)  is  particularly  useful  and  suggestive. 


3.  A  Conservation  Law  with  Fading  Memory.  We  study  the  model  nonlinear,  history-boundary 
value  problem 

t 

(3.1)  u*.  +  w(u)  +  /  a'  (t  -  T)*{u(T,x>>  at  -  f(t,x)  <-«  <  t  <  -,  0  <  x  <  1)  , 

_®  x 

subject  to  the  periodic  boundary  condition 

(3.2)  u(t,0)  -  u( t , 1 )  ) 

the  history  of  the  solution  u  is  assumed  to  be  known  up  to  time  t  »  0,  i.e. 

(3.3)  u(t,x)  -  v(t,x)  (-«•  <  t  <  0,  0  <  x  <  1)  , 

where  v  is  a  given  smooth  ( c 1  <  ( — 0  ]  x  (0,1]))  function  which  satisfies  (3.1),  (3.2) 
for  t  <  0.  In  (3.1)  :  R  ♦  R,  f  :  R  *  [0,1]  ♦  H,  and  a  s  [0,")  ♦  R  are  given 

functions  satisfying  assumptions  analogous  to  those  for  the  viscoelastic  problem  outlined 
in  Section  2.  Our  task  in  this  section  is  to  determine  a  smooth  extension  u(t,x)  of 
v(t,x)  on  (-“,")  X  [0,1]  which  satisfies  (3.1)  and  (3.2),  and  to  study  the  asymptotic 
properties  of  u  as  t  ♦  “.  In  order  to  do  this  the  history  v  and  the  forcing 
function  f  will  have  to  be  taken  sufficiently  smooth  and  "small"  in  suitable  norms. 

The  requirement  that  the  history  function  v  should  satisfy  (3.1)  and  (3.2)  for 
t  <  0  is  motivated  by  the  viscoelastic  problem  described  in  Section  2.  The  history  value 
problem  in  which  the  function  v  is  sufficiently  smooth  for  t  <  0  (but  need  not  satisfy 
(3.1))  is  also  of  interest,  and  can  be  studied  by  the  same  methods  (see  further  remarks 
below) . 

The  basic  assumptions  for  the  global  existence  theory  are  as  follows.  Concerning  the 
constitutive  functions 

(c)  w,i|>  e  c2<R),  w(0)  =  W0)  =  o,  w'(0)  >  o,  <t”(0>  >  o  , 

concerning  the  kernel  a: 

2  1 

(a)  a  e  W  '  (0,°*)  and  a  is  strongly  positive  on  [O,09)  ; 

concerning  the  forcing  term  f: 


f»ft,fx  e  c[(-“,«)il2(o,1))  n  l2k--,«) »l2(0,i>]  , 
f(t,x)  »  f^t.x)  +  f2<t,x)  , 

(£)  ■  f1tt,f2tx  e  l2[(-«,-)il2(0,di  , 

1 

/  f(t,x)dx  -  0  {-*•  <  t  <  «)  , 

0 

and  to  measure  the  "size"  of  the  forcing  term  we  define 

1  •  1 

F  -  sup  /  {f2  +  f2  +  f2}(t,x)dx  +  /  /  {f2  +  f2  ♦  f2  ♦  f2  +  f2tx}dxdt  » 

(--,»)  o  -“0  c 

concerning  the  history  vs 

v,vt,vx,vtt'vtx,vxx  e  c((-"»“)>l2(0,1)]  n  L2 [ »L2(0, 1 )] , 

1 

(H)  /  v(t,x)dx  -  0  (-«•  <  t  <  0)  , 

0 

and  v  satisfies  (3.1),  (3.2)  for  t  <  0  . 

Analogous  to  the  "equilibrium  stress"  for  the  viscoelastic  body  we  define  the  constitutive 
function  x  :  R  +  R  by 

XC)  =  v5(«)  -  a(0)t|>(  ’ )  , 

and  we  assume  that 

(3.4)  X’(0>  -V’(0)  -a(0)*'(0)  >  0  . 

The  following  equation,  obtained  from  (3.1)  by  carrying  out  an  integration  by  parts  with 
respect  to  T,  will  play  a  crucial  role  in  the  analysis: 

t 

(3.5)  ut  +  X(u)  +  /  a(t  -  T)!|>(u(T,x))  _dT  =  f(t,x)  (-“  <  t  <  »,  0  <  x  <  1)  . 

x  _<•  xT 

It  is  clear  that  the  problem  consisting  of  (3.5),  (3.2),  (3.3)  is  equivalent  to  the 
original  problem  (3.1),  (3,2),  (3,3). 
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Upon  setting 


0 

(3.6)  h(t,x)  -  -  /  a*  (t  -  XJWvlT.xn^dx  +  f(t,x)  (0  <  t  <  “,  0  <  x  <  1) 


(3.7)  Ug(x)  -  v(0,x)  (0  <  x  <  1)  , 

the  history-boundary  value  problem  (3.1)— (3.3)  reduces  to  the  lntia 1-boundary  value  problem 


(3.8) 


t 

“t  *  ^(u,x  *  /  »'(t  •  T)*(u(T#x))  dt  -  h(t,x)  (0  <  t  <  “,  0  <  x  <  1) 

0  * 
u(t,0)  »  u(t,1)  (0  <  t  <  “) 

u(0,x)  -  Uq( x)  (0  <  x  <  1)  , 


1 

where  /  h(t,x)dx  «  0  (0  <  t  <  “)  [use  v(t,0)  -  v(t,1),  and  assumptions  (f)]. 

0  1 

Conversely  a  solution  of  (3.8),  where  /  uQ(x)dx  “  0,  can  be  shown  to  solve  (3.1)-(3.3) 

0 

by  constructing  a  smooth  function  v(t,x)  on  <-“,01  *  [0,1],  satisfying  (3.2), 


1 

/  v(  t ,  x )  dx  «  0 
0 


(3.9) 


(-“  <  t  <  0),  and  also  requiring  v(0,x)  »  uQ(x),  as  well  as 
vt(0,x)  «  -v(ufl(x))x  +  h(  0,x)  -  u^x) 

vtt(0,x)  -  -v"(u0(x))uj(x)u1(x)  -  v>'(u0(x))uj(x) 

-  a'lOJWUglx))^  +  hfc(0,x)  , 


and  defining 


(3.10)  f(t,x)  - 


*t  +  *(v)x  +  /  a'(t  -  T)<i(v(T,x))  dx  (-“  <  t  <  0,  0  <  x  <  1) 


h(t,x)  +  /  a'  ( t  -  T)i|»(v<x,x) )  dx  <0  <  t  <  -,  0  <  x  <  1)  ; 


the  requirements  (3.9)  insure  that  f  defined  by  (3.10)  has  the  necessary  smoothness 
properties  across  t  «  0  required  In  the  existence  theory. 

The  main  global  result  Is 

Theorem  3.1.  Let  the  assumptions  (a),  (c),  and  (3.4)  be  satisfied.  There  exists  a 
constant  U  >  0  such  that  for  every  forcing  term  f  satisfying  assumptions  (f)  with 


2 

F  <  M  ,  and  for  any  history  function  v  on  (-",0}  *  [0,1]  satisfying  assumption*  (H), 
there  exists  a  unique  function  u  on.  (-*,*)  *  [0,1]  with 

u'ut'ux,utt'utx'uxx  e  c[ (■"'“)  »L2(0,D]  n  L2[ (-»,«>  >L2< 0,1)]  satisfying  <3.1)-(3.3>,  as 
1 

well  as  /  u(t,x)dx  -  0  ( <  t  <  ••) .  Moreover, 

0 


(3.11)  u(t,x),  ut(t,x),  ax(t,x)  *0  as  t  ♦  •  , 

uniformly  on  [0,1]. 

The  proof  of  Theorem  3. 1  will  be  given  In  Section  4.  It  uses  the  general  strategy 
developed  In  [5] ,  although  there  are  technical  differences  In  details.  One  first 
establishes  the  existence  of  a  local  solution  u  on  a  maximal  Interval  (-*',Tq)  *  [0,1], 
with  the  property  that  when  Tg  <  "  a  certain  norm  of  u  becomes  infinite  as  t  +  Tg 
(see  Proposition  4.1  below).  One  then  uses  a  combination  of  energy  methods  and  properties 
of  strongly  positive  kernels  to  show  that  the  integral  in  (3.1)  exerts  a  dissipative  effect 
resulting  in  the  aforementioned  norm  remaining  uniformly  bounded,  independent  of  Tg, 
provided  the  constant  li  in  Theorem  3.1  is  sufficiently  small.  ThUB,  in  particular 
Tg  »  +”  and  the  smooth  solution  exists  globally. 

Remark  3.2.  It  follows  by  standard  regularity  techniques  that  the  solution  u  of  (3.1), 
(3.2),  (3.3)  is  C1  smooth  on  R  x  [0,1], 

Remark  3.3.  A  result  similar  to  Theorem  3.1  (established  by  the  same  methods)  holds  for 
the  boundary-initial  value  problem 


(3.12) 


u..  +  ^(u)„  +  /  a'(t  -  x)  ij)(u(  x,x) )  dx  »  h(t,x)  (0  <  t  <  *,  0  <  x  <  1)  , 
0 

u( t ,0 )  -  u ( t , 1 )  (0  <  t  <  1)  , 

u(  0 , x)  “  uo<*)  (0  <  x  <  1,  Ug(0)  «*  u g  ( 1 ) )  > 


here  v“,  and  a  satisfy  the  assumptions  of  Theorem  3.1,  h  defined  on  [0,®)  x  [0,11 


satisfies  assumptions  (f)  modified  in  the  obvious  way,  and  the  initial  datum  uQ  e  HMO,  1), 
1 

with  /  u0(x)dx  "  0.  For  the  theorem  to  hold  one  must  require  that  lu  I  ,  and  a 
0 

suitable  norm  of  h  are  sufficiently  small,  as  can  be  seen  from  a  detailed  examination  of 


the  corresponding  estimates.  (See  Remark  4.3  below) 


Remark  3 .4 .  The  requirement  that  the  history  v  satisfies  (3.1),  (3.2)  for  t  <  0,  and 

2 

the  condition:  the  constant  U  (where  F  <  U  )  of  Theorem  3.1  is  "small",  imply  that  the 
history  function  v,  as  well  as  the  forcing  term  f  in  (3.1),  are  both  small  in  suitable 
norms . 

Remark  3.5.  If  \f  =  in  (3.1),  Theorem  3.1  can  be  applied  without  any  change,  provided 
0  <  a(0)  <  1  (in  order  that  (3.4)  is  satisfied).  However,  in  this  special  case  the 
somewhat  different  energy  techniques  of  Dafermos  and  Nohel  [4] ,  or  of  Staffans  [20]  ,  or  the 
method  of  Riemann  invariants  of  MacCamy  [11]  can  also  be  used. 

Remark  3.6.  It  will  follow  from  the  proof  of  Theorem  3.1  in  Section  4  that  because  the 

problem  (3.1)-(3.3)  is  in  one  space  dimension  and  on  a  finite  space  interval  it  is 

sufficient  to  obtain  global  estimates  for  the  derivatives  (in  this  case  estimates  of  utt> 

•  2  2  2 

utx'  “xx'  b“cause  (3.1)  is  of  first  order)  in  the  L  (L  (0,1))  and  L  (L  (0,1))  norms, 
and  then  to  make  use  of  the  Poincare  inequality  to  estimate  lower  order  terms.  However, 
one  cannot  apply  this  method  to  obtain  global  estimates  for  the  pure  Cauchy  problem 
consisting  of  the  Volterra  equation 

t 

(3.13)  u..  ♦  *(x)  +  /  a’(t  -  t)*(u(t,x)>  dr  -  h(t,x>  (0  <  t  <  -,  x  e  R)  , 

0  * 

together  the  initial  condition  u(0,x)  -  uQ(x),  x  &  R,  because  x  e  R  and 
the  Poincare  inequality  does  not  apply.  For  the  local  existence  result  (analogue  of 
Proposition  4.1)  one  can  circumvent  this  difficulty.  The  same  oomments  apply  to  the 
analogous  pure  Cauchy  problem  associated  with  the  second-order  nonlinear  Volterra  equation 
(2.8).  If  p  =  p  in  (3.13)  or  (2.8)  the  analogous  pure  Cauchy  problems  can  be  treated  by 
either  the  methods  of  [4],  [11]  or  of  [20],  because  global  estimates  of  u  and  of  the 
derivatives  of  u  in  appropriate  norms  are  obtained  successively  from  the  equations . 
However,  if  ♦  J  <)i  in  (3.13)  or  the  pure  Cauchy  problem  associated  with  (2.8)  these 
mathematically  interesting  problems  remain  to  be  tackled. 
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4.  Proof  of  Theorem  3.1)  a.  tool  Theory.  This  is  carried  for  the  houndary-initial  value 
problem  (which  was  shown  to  be  equivalent  to  the  history  value  problem  (3.1 )  —  (3.4)  in 
Section  3  for  a  suitable  choice  of  the  forcing  term  f  -  see  (3.10))t 


u,.  +  v>(u)  +  /  a'(t  -  t)(i(u(t,x))  dt  »  h(t,x)  (0  <  t  <  ",  0  <  it  <  1) 

0  x 
u(t,0)  »  u(t,1),  0  <  t  <  •  » 

u(0,x)  -  un(x),  x  e  [0,1]  . 


He  make  the  following  assumptions! 


'.a-  e  l’io,-),  e  c2 (*) ,  *(0)  -  *(0)  -  o,  *'(o)  >  o  , 


there  exists  a  constant  <  >  0  such  that  ^ 1 ( 5 )  *  *  >  0  (5  @  R) , 


h  s  [0,«]  *  [0,1]  ♦  R  ,  h( t ,x)  -  h,(t,x)  +  h2(t,x),  h,ht,hx  B  C([0,-)»L  (0,1)), 

/  h(t,x)dx  -  0  (t  e  (0,-)),  n1tt,h2t  6  L2([0,-)»L2(0,D),  and  u„«  HJ(0,1), 

01 

/  ug(x)dx  -  0.  The  reader  will  note  that  the  full  strength  of  assumptions  (a)  and 
0 

assumption  (3.4)  are  not  used  in  Proposition  4.1  below;  on  the  other  hand  <f> '  is  now 

required  to  be  bounded  away  from  zero  (oompare  with  assumptions  (c))« 

Proposition  4.1.  tinder  the  above  assumptions  there  exists  a  0  <  Tg  <  "  and  a  unique 

function  u  e  c'dO.Tg)  *  [0,1]  with  utt,utx,uxx  e  C(  [0  ,T]  »L2(0, 1 ) )  for  every 

0  <  T  <  Tg,  such  that  u  satisfies  (4.1)  on  [0,Tg)  x  [0,1]  and  /  u(t,x)dx  -  0. 

0 

Moreover,  if  Tg  <  • 


(4.2)  sup  /  [u2(t,x)  +  u|(t,x)  +  •••  +  u2x(t,x)]dx  -  . 

«;• 

It  is  clear  that  with  h  defined  (3.6),  uQ  by  (3.7)  the  problem  (3.1)— (3.3) 
satisfies  the  assumptions  of  Proposition  4,1. 

Remark  4 .2  .  A  similar  result  holds  for  the  pure  Cauchy  problem  associated  with  (4.1),  i.e. 
no  boundary  condition  is  specified  and  u(0,x)  -  Ug(x),  uQ  B  H  (R);  however,  the  function 
space  X (M,T)  below  must  be  specified  differently. 

The  proof  of  Proposition  4.1  is  very  similar  to  that  of  Theorem  2.1  in  [5],  and  will 


only  be  sketched.  Let  M  >  0,  T  >  0  and  let  X(M,T)  denote  the  set  of  functions 
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w(t,x)  on  [0,T]  *  [0,1]  with  w,wt.wx,wtt,wtx,wxx  e  C( (0,T;L2(0, 1 ) ) ,  w(t,0>  -  w(t,1), 

w(0,x)  «  Ug(x),  /  w(t,x)dx  »  0,  /  w.(t,x)dx  «  0,  0  <  t  <  T,  and 
0  0 

1 

(4.3)  sup  /  [w£t(t,x)  +  w^x(t,x)  +■  wxx(t,x)]dx  <  M2  . 

[0,T]  0 


It  follows  from  the  Poincare-type  inequalities  (see  application  of  Lenmas  A.1  and  A.2  in 
Appendix)  and  from  (4.3)  that  if  we  X(M,T)  then 

(4.4)  w2(t,x)  +  w2(t,x)  +  w2(t,x)  4m2  (0  <  t  <  T,  0  <  x  <  1)  . 

Let  S  :  x (M,T)  *  C^dO.T]  x  [0,1]  be  the  mapping  which  carries  w  e  X(M,T)  into  the 
solution  of  the  linear  problem  consisting  of 


t 

(4.5)  u  +  *>’(w)uv  «  -  /  a*  <t  -  T)*(w(T,x))  dt  +  h(t,x)  (0  <  t  <  T,  0  <  x  <  1 )  , 
6  *  *  X 


and  of  the  boundary  and  initial  conditions  in  (4.1).  It  is  clear  that  a  fixed  point  of  S 

is  a  solution  of  (4.1)  on  [0,T]  x  [0,1].  Also  note  *'(w)  is  W1'  smooth  and  ¥>'(w)t 
■  2 

and  *'(w)x  are  in  L  ([0,T];L  (0,1)).  Moreover,  if  g(t,x)  denotes  the  right-hand  side 

of  (4.5),  g(t,x)  ■  g^ (t,x)  +  g2(t,x),  then  g  satisfies  the  same  assumptions  as  h 

does,  and  by  fairly  standard  theory  for  linear  problems  utt»utx,uxx  e  C{ [0,T] ;L2(0, 1 ) ) i 

embedding  type  arguments  then  yield  that  u  e  c'dO.T]  x  [0,1]).  Thus  it  suffices  to  show 

that  the  map  S  has  a  unique  fixed  point  u  in  X(M,T)  .  Once  this  has  been  demonstrated 

1 

it  follows  from  the  assumption  /  h(t,x)dx  »  0,  t  8  [0,1],  the  equation  in  (4.1),  and 

0 

3  t  1 

from  the  boundary  condition  that  *-  /  u(t,x)dx  *>  0  for  t  e  [0,TJ»  since  /  uQ(x)dx  »  0, 
1  *  0  0 
one  then  also  has  /  u(t,x)dx  «  0,  t  e  [0,T]  . 

0 

The  remainder  of  the  proof  of  Proposition  4.1  is  completed  in  the  following  steps: 

(1)  analogous  to  the  proof  of  Lemma  2.1  in  [5]  (here  the  argument  is  shorter),  use  the 
standard  energy  method  for  linear  problems  to  show  that  when  M  is  sufficiently  large 
and  T  >  0  is  sufficiently  small,  S  maps  X(M,T)  into  itself:  (li)  define  the  metric 
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1 

P(u,u)  -  max  {/  t(ut  -  ut)2  +  (u  -  u„)2] (t,x)dx  , 

[0,T]  O' 

where  u,u  e  X(M,T);  by  the  lower  semloontinuity  of  norma  in  a  Banach  space,  X(H,T)  ia 
complete  under  Pj  analogous  to  Lemma  2.2  In  [5]  show  that  S  is  a  strict  contraction 
on  X(M,T)  under  the  metric  P;  (iii)  by  Banach's  fixed  point  theorem  the  map  S  has  a 
unique  fixed  point  u  e  C^([0,T]  *  [0,1])  for  M  sufficiently  large  and  T  sufficiently 
small,  which  solves  (4.1);  the  existence  of  the  maximal  interval  of  validity 
[0 ,Tq )  «  [0,1]  of  the  solution  u  satisfying  (4.2)  is  established  in  a  standard  manner 
(see  [5]).  This  completes  the  sketch  of  the  proof  of  Proposition  4.1. 

_b .  Global  Theory.  By  the  constitutive  assumptions  (c)  and  (3.4)  there  exist 
5  >  0 ,  x  >  o  such  that 

(4.6)  *'<5)  >  x,  *'(5)  >  x,  x'  (£)  >  <  (151  <  4)  . 

If  necessary  modify  outside  [-4,4]  such  that  s 9  6  C2(*)  and  f'(  5)  >  x  (5  e  R)  .  To 
prove  the  existence  of  the  global  solution  u  on  (-**,0*)  x  [0,1]  of  the  history  value 
problem  (3.1)-(3.3)  asserted  in  Theorem  3.1,  let  u  be  the  unique  solution  on  the  maximal 
interval  (*-,,T(J)  x  [0,1]  guaranteed  by  Proposition  4.1,  and  assume  that  0  <  T0  <  “•. 

For  0  <  T  <  T„  let 

1 

(4.7)  U(T)  -  sup  /  [u2(t,x)  +  u^(t,x)  +  u2(t,x)  +  •••  +u2x(t,x)]dx 

t--,T]  0 

T  1 

+  /  /  [u2  +  u?  +  u2  +  •••  +  u2  1  dxdt  , 

-“0 

where  •••  stand  for  the  terms  u£t  and  u£x .  Recall  that  TQ  is  characterized  by  (4.2) 
and  thus  the  first  Integral  in  U(T)  tends  to  infinity  as  T  ♦  TjJ;  also  recall  that 
u  =  v  for  t  <  0.  The  basic  strategy  of  the  proof  is  the  same  as  in  [5]  >  we  wish  to  show 
that  there  exist  constants  v  >  0  (v  <  4  in  (4.6))  and  X  >  0,  independent  of  T,  such 
that  if 

(4.8)  I u ( t , x) | 2  +  |ux(t,x)|2  ♦  |ut(t,x)|2  <  v2,  v  <  S  , 
on  <-“,T]  *  [0,1]  then 
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•  wu 


I 


if.; 

* 


(4.9)  U(T>  <  KF  , 

where  F  is  the  constant  defined  under  assumptions  (f).  The  proof  that  (4.8)  implies 

(4.9)  will  be  outlined  below  using  energy  estimates.  Once  this  claim  is  established  the 

proof  of  Theorem  3.1  is  completed  as  follows.  First,  by  the  assumptions  on  the  history 
v(t,x),  (4.8)  holds  as  a  strict  inequality  for  t  >  0  sufficiently  small.  Next,  by 

the  Poincare  inequality  (see  Lemma  A. 1  in  Appendix)  and  the  definition  of  U(T) 

(4.10)  |u(t,x)|2+  |ux(t,x)  |2  |ut(t,x)|2 

<  /  [u2(t,x)  u2  (t,x)  +  u?  (t,x))dx  <  U(T) 


2  V  2 

on  (-*»,T]  *  [0,1].  Choosing  the  constant  V  <  —  and  F  <  u  ,  (4.10)  shows  that  (4.9) 

K  2 

implies  (4.8)  (as  a  strict  inequality).  Therefore  if  F  <  u2  <  the  estimates  (4.8) 

and  (4.9)  both  hold  for  every  T  8  (-“,TQ),  and  consequently  by  Proposition  4.1  (see 

especially  (4.2))  TQ  -  +“.  Moreover,  (4.7)  and  (4.9)  imply  that 

(4.11)  u'ut'ux«utt'utx'uxx  8  t*((--,*)>L2{°,1))  H  L2( (-»,») iL2<0.1))>  . 

2  2 

But  u,ut,ux,utt, ufcx , uxx  e  L  ( ( -"*,*•)  tit  (0,1))  also  implies  that 


(4.12) 


u,ut,ux  *0  in  L*(0,1)  as  t 


which  in  view  of  u,ut,...,u  @  b  ((-*V*)>b  (0,1))  yields  (3.11)  and  completes  the  proof. 

It  remains  to  establish  that  (4.8)  implies  (4.9).  For  this  purpose  we  will  need  the 
following  properties  of  strongly  positive  Xernels.  Introduce  the  notation 

t 

(a*g) (t)  -  /  a(t  -  T)g(T)dx 


and 


Q[a,w>s]  »  /  w(t) (a*w) (t)dt  . 


Let  assumptions  (a)  of  Theorem  3.1  be  satisfied.  There  exist  constants  B,  Y  >  0  such 

2 

that  for  every  s  @  *  and  for  every  w  e  L  t-*,s) 
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xr 


(4.13) 


/  [(a*w)(t)l2dt  <  BQ[a,wis] 


B  -  \  'a'2. 


2i  +i'*,,ai 

l'(0,-)  “  L  (0,«) 


(4.14) 


/  [(a'*w)(t)]2dt  <  YBIa,w»s]  , 


Y-S  u*,21 

T.  I  ft  .«•! 


2,  -  i  '*"«2 

L  <0,«)  L 


1  * 

(0,“) 


The  eat lm tea  (4.13),  (4.14)  which  have  alao  been  used  in  [5]  are  essentially  contained  in 
Staffans  [20].  Another  important  property  deals  with  the  resolvent  kernel  k  of  the 


linear  Volterra  operator 


•fi' <0)y  +  ♦' (0)a'»y 


defined  to  be  the  solution  of  the  linear  Volterra  equation 


*'(0)k(t)  +  *' (0)  /  a'(t  -  T)k(T)dT  -  -+• (O)a' (t)  0  <  t  <  •  . 

0 


Lemma  4.2.  If  assumptions  (a)  are  satisfied,  and  if  ^'(0)  >  0,  <>'(0)  >  o,  then 

k  e  l’(0,-). 

The  proof  of  Lemma  4.2  is  given  in  [51. 

The  first  estimate  needed  for  the  proof  that  (4.B)  implies  (4.9)  is  obtained  from 
equation  (3.5).  [Recall  that  (3.5),  (3.2),  (3.3)  Is  equivalent  to  the  original  problem 
(3. 1)— (3.3) ,  and  it  is  assumption  (3.4)  concerning  x  which  will  play  an  important 
role.]  Multiply  (3.5)  by  i>(u)  •nd  Integrate  over  t-“,s]  *  [0,1],  s  <  TQ.  After 

several  integrations  by  parts  in  which  the  boundary  condition  is  invoked  we  obtain 
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(4.15) 


4  /  X'  <u(s,x)  M>'  (u(a,x))u*(a,x)dx  +  Qta.s^fu)  ) 

*  0  *  xt 

.  3  1 

•  —  •?  S  I  (X1  (uH>"  (u)  +  X“ Cu) (u)]u  u  dxdt 
--  0  r 

1  8  1  1 

-  r  /  /  **(u)u  u^dxdt  +  /  f  (8,x)<Mu(s,x)  )  dx 

-•  0  e  0  * 

8  1 

-  /  /  f.  ♦(u)  dxdt  . 

‘  '  t  x 

0 

In  oontraat  to  the  analogous  calculation  in  [5,  see  (3.21)),  no  useful  information  is 
extracted  here  from  integration  of  the  term  xt  over  (-“.si  *  [0,1], 

Remark  4.3.  When  obtaining  the  analogous  estimates  for  the  boundary-initial  value  problem 
(3.12)  (see  Remark  3.3),  the  analogue  of  equation  (3.5)  from  which  the  global  estimates  are 
calculated  is 


t 

ut  +  X(u)x  +  /  a(t  -  T)f(u(T,x)  )xTdT  -  h(t,x)  -  a(t)i|i(u0)  . 

0  x 


To  simplify  several  of  the  estimates  which  follow  we  make  the  additional  assumption 

that  ♦  (and  hence  also  x)  @  C^(*)j  the  alternative  Is  to  employ  difference  quotients 

and  pass  to  limits  as  in  (5).  Differentiate  (3.5)  with  respect  to  t  (use 

(a*g)(t)  “  /  a(S)g(t  -  S)d£,  differentiate,  and  then  change  variables)  obtaining 
0 

(4.16)  utt  +  Xlu)^  +  a**(u)xtt  -  ffc  (_«  <  t  <  «,  0  <  x  <  1)  . 

Multiply  (4.16)  by  an<5  integrate  over  (-*,s)  *  (0,1).  After  several 

integrations  by  parts  the  result  of  this  tedious  calculation  is 
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(4.17) 


a  1 


Ja  -  /  /  X '  ( u)  *•  ( u) u^u^dxdt 


^  /  X'<u(a,x))i(i,(uf»»x>)u^t<s,x)dx 


a  1 


~  J  f  /  lx"(un|i,(u)  +  X* < «)♦"(«> 3 u  u  dxdt 
-»  n  xt  c 


which  respectively  give  the  first  terra  on  the  left  aide  and  last  tern  on  the  right  side  of 
(4.17), 


J1  "  /  /  X"(u)4i"'(u)u^u^dxdt  , 

-m  0  x  c 

si 

•*2  “  2  /  /  X“(u)  ♦"(ulu  u^u  dxdt  , 

-»  0  xc  r  x 

ra  r1  2 

J3  *  /  /  X"(u) V"(u)uttuxufcdxdt  , 

1 

J4  -  /  X"(u(s,x) )♦' (u(s,x) )u  <s,x)u  (s,x)u  (s,x)dx 

0  Z  X  xt 

si  si 

-  /  /  x"<«)+'(u)u  u  dxdt  -  /  /  u  u  [x"(u)**(tt)u  ]  dxdt 

_«•  o  *>-  t  _»  q  xt  x  t 


s  1 

J5  ■  /  /  X'(u)i|r" '(uju^u^u^dxdt  , 


0 

s  1 

/ 
0 


J6  *  2  /  /  X'(u)'l»"(u)ux  u  dxdt  , 

0 


s  1 


J7  -  /  /  X’(u)*"(u)u  ,u..u  dxdt  . 

7  0  xt  tt  x 


The  next  estimate  follows  from  the  identity 

(4.18)  a(0)<|>(u)  »  a'*<((u)  +  a*ip(u)  . 

xt  xt  xtt 

which  is  derived  by  integrating  a*'(i(u)  by  parts.  Multiply  (4.18)  by  uvf.  and 

Xtt 
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I 


-> 


I  , 


r  1 

i  ■ 


integrate  over  (-**,s)  *  (0,1)  »  then  in  each  term  on  the  right  hand  side  use  the  Cauchy 
Schwartz  inequality  and  (4.13),  (4.14).  This  leads  to  the  estimate 


s  1 

(4.19)  a(0)  /  /  u  t(u)  .dxdt 

-«  o  tx  ** 


<  If  f  -^dxdt]  ^  {(rOta,s»Wu)v.)>  1/2  +  (BQ[a,a,*(u>  .  .]>  '/2  } 
-»  0 


xt' 


'  xtt J 


Using  «(u)xt  -  (u)uxt  +  ♦"(uju^  and  (4.6)  on  the  left  aide  of  (4.19)  and 
a2  C  2 

eb<2^  +  '2b,e  -  2/«a(0)  on  the  right  aide  of  (4.19)  givea  the  useful  estimate 


(4.20) 


u*xdxdt  -  Qta.sjD.tu)^] 


B  8  1 

icaToT  Ota, ai <>(0)^1  <  a(0>  /  /  \|>"(u)u  u  u  dxdt 

-•0 


Sext  write  (4.16)  in  the  form 

utt  -  -<X"(u)uxut  +  xMulu^  +  +  ft  i 


square  both  sides,  integrate  over  (-“.s]  *  (0,1), 
use  (4.13)  to  obtain  the  estimate 


use 


(  I  a.)2  <  n  l  a2  (n 
i-1  i=*1 


4 ) ,  and 


si  si 

(4.21)  /  /  u£tdxdt  -  4  /  /  (X'(u)]2u2tdxdt 


•  0 


-«  0 


8  1 

-  48Q[a,si*(u)  1  «  4  /  /  (x"(u)]2u2u2dxdt  +  4  /  /  f^dxdt 

-»  o  --  0 


2  2  2. 


8  1 

/ 
o 


Me  now  return  to  (3.1),  differentiate  with  respect  to  t  and  obtain 

“tt  *  -<*’(u>«xt  +V>"(u)uxut  ♦  a*«»Mu)Uxt  +  a***-(u)uxut)  +  ffc  , 
square  both  sides,  integrate  with  respect  to  x  over  [0,1],  and  evaluate  at  t  »  s. 
This  gives  the  estimate 
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I 


■f 


i  i 

(4.22)  /  u£t(s,x)dx  -  5  /  (*>'(u(s,x))]2u2x(s,x)dx 

0  0 


-  5  la' i  sup  /  ($' (u(t,x) ) !  u  ( t ,  x ) dx 

l’<0,-)  (--,8]  0  tx 

1  1 
<  ?  /  {^"(u(s,x) ) j 2u2( s,x)u?(  s,x) dx  +5  /  f2(s,x)dx 

0  0 


♦  5  la* I2  sup  /  («"(u(t,x)J 2u2(t,x)u^(t,x)dx 

L’CO,-)  -«<t<8  0  * 


Next,  differentiate  (3.1)  with  respect  to  x  obtaining 

utx  +  *<u>xx  +  a'**(u)xx  ’  fx  ' 

and  write  ^  ( u>Xx  ”  ^'*u*uxx  +  “  ^'(0)uxx  +  ~  v>'(0)]uxx  +  ^"fu^ux»  and 

similarly  for  ^u^xx*  This  gives  the  equation 

(4.23)  if'  (0)u„„  +  ip'  (0)a'*u  -  X(t,x)  , 

**  XX 


X(t,x)  »  -utx  +  fx  -  (^'(u)  -  *  C  0 )  I  v»xJt  -  v>"(u)ux 

-  a'*[ip'(u)  -  *'<0)]u  -  a'*rtu)u^  . 

Letting  k  he  the  resolvent  of  the  operator  on  the  left  side  of  (4.23)  we  have 


(4.24) 


^•(0)uxx(t,x)  -  X(t,x>  +  (k*X) ( t,x)  . 


By  Lemma  4.2  k  e  L  (0,"),  and  this  gives  the  estimates 


(*>'(0)]2  /  u2x(s,x)dx  <2  /  X2  ( s,x)  dx  +  2  »k  I2  sup  /  X2(t,x)dx  , 
0  0  L  1  (  1 8 }  0 


31  s  1 

(V>'(0)]2  /  /  u2xdxdt  <  2  [  1  +  Ik  I  .1  /  /  X2dxdt 

-»  0  l'  -•  0 


To  complete  the  proof  that  (4.9)  implies  (4.9),  we  first  use  (4.8),  (4.6),  (f),  and 
the  fact  that  for  |u(t,x)|  <  v,  the  derivatives  1*0  '  (u(  t,x) )  I ,  |sj"(u(  t,x) )  |  , 


1 


t 


1 


k',,(u(t,x))|,  |*,<u(t,x))|,  |*"(u(t,x))|,  |*,"(u(t,x))|,  lx*  (u(t,x) )  |  /  Ix"(u(t,x))|, 

I X' * ' (u(t,x) ) |  are  bounded  by  a  constant  C  >  0  In  order  to  simplify  the  basic  estimates 
derived  above;  here  we  have  used  the  additional  simplifying  assumption  that  and  hence 

also  X  e  c3(R);  this  can  be  avoided  as  in  [5]. 

Thus  (4.15)  becomes 


7-  /  ux(s ,x) dx  +  Q(a,*(u)  »s)  <  §•  vc 2  /  /  u2dxdt 

4  0  xz  i  --  0 

+  /  /  u|dxdt  +  Fc[^r  +  — r~r  ]  ; 

2  0  **  2C  /v 


similarly  (4.17)  simplifies  to 


<2  f1  2 

—  /  u2t(s,x)dx  +  Qla,1>(u)xtt»s) 


1  1 
<  c2v  /  u2t(s,x) dx  +  cv[-l  +  c]  /  u2(s,x)dx 


si  si 

CV[F  +  C  +  CV  +  CV2]  /  /  u?dxdt  +  c2v2  /  /  u?dxdt 

-«  0  --  0 


+  cv[j  +  (f  +  |  vjc]  /  /  u2xdxdt 

-•0 


cv[2  +  |  c  +  (|  +  c)v]  /  /  u2tdxdt  +  cf(v  ♦  ~  + 


the  estimate  (4.20)  becomes 


1  >  r.**-* -•*■»**•*  *-vr*v»-4 


I 


4 


(4.29)  /  u^dxdt  -  Qla.s^u)^] 


B  I  B  I  . 

Qta.af^u)^!  <  a(0)  (/  /  u^dxdt  +  /  /  lylxdt); 


estimate  (4.21)  simplifies  to 


(4.30) 


si  si 

/  /  u£tdxdt  -  4C 2  /  /  u£fcdxdt 

-m  0  -*•  0 


-  40Q[a,s>'t>(u)  1  <  4CV  /  /  u£dxdt  +  4F  » 

xrt  -•  0 


estimate  (4.22)  simplifies  to 


(4.31) 


/  u|t(s,x)dx  -  SC2  /  u2Jt(s.x)dx 
0  0 

1  2  2  1 

-  5C2la’l  .  sup  /  U2  (t,x)dx  <  5CTv2  /  u2(s,x)dx 
L1  -«<t<8  0  0 

1 

+  5C2v2 la* •  sup  /  u2(t,x)dx  ♦  5F  > 

L1  -«*<t<s  0 


the  estimate  (4.25)  becomes 


(4.32) 


'2  /  u2  (s,x)dx  <2  /  X2(s,x)dx  +  2lkl2  sup  /  X2(t,x)dx  ; 

•  -  *  ,  —  .1  St 


L  (-“,S]  0 


finally,  (4.26)  becomes 


SI  *  1  * 

/  /  u2x<Jxdt  <  2{1  *  Ikl  j  /  /  X2dxdt  , 

-»  o  l  -•  n 


where  in  the  estimates  (4.32)  and  (4.33)  X(t,x>  is  the  function  (depending  on  n,  ux, 
utx,  a,  and  f)  given  preceding  formula  (4.24). 
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We  now  focus  our  attantlon  on  tha  simplified  estimates  (4.27)  through  (4.33). 

By  tha  Poincare  inequality  (see  Appendix)  o(T) ,  given  by  (4.7),  can  be  majorised  by 


(4.34)  U(T)  <  sup  [/  (2ux(t,x)  +  u|t(t,x)  ♦  2u?  ( t,x)  +  u  (t,x)]  ]dx 

(--,T)  0  XX 

T  1 

+  /  /  (u?t  ♦  2u?  ♦  3uJx]dxdt  . 

-»  0 

Moreover,  each  term  on  the  right  hand  aide  of  (4.34)  can  be  majorised  by  a  suitable  linear 
combination  of  the  left  hand  aides  of  tha  estimates  (4.27-4.33).  On  the  other  hand,  each 
term  on  the  right  hand  sides  of  the  estimates  (4.27-4.33)  can  be  majorised  by  terms  of  the 
form  0(F),  or  0(v)u<T),  or  eu(T)  ♦  c(C)F(T)  for  any  e  >  0»  the  last  of  these  comes 
from  estimating  the  right  hand  sides  of  (4.32)  and  (4.33).  This  combined  with  (4.34) 
yields  the  final  estimate 

(4.35)  U(T)  <  <0( V)  ♦  0(e>}0(T)  +  c(C)P  . 

Therefore,  fixing  v  >  o,  c  >  0  sufficiently  small,  (4.35)  yields  (4.9),  assuming  that 
(4.8)  holds.  This  completes  tha  proof  of  Theorem  3.1. 
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5.  Development  of  Singularities 


Ha  consider  the  pure  Cauchy  problem 


u,.  +  tf(u)  +  a'**(u)  -  0 

(5.1)  {  x 

u(0,x)  -  Ug(x)  (x  e  R)  > 

throughout  this  section  *  will  denote  the  convolution  on  [0,t]  (not  (-*,t] ) .  For  a 

discussion  of  the  existence  of  smooth  solutions  of  (5.1)  on  (0.**)  *  R  for  sufficiently 

smooth  and  ’small"  data  we  refer  the  reader  to  Remarks  3.3,  3.6.  It  is  known  (see 

Proposition  4.1  and  Remark  4.2)  that  if  a,  <t,  t  satisfy  the  assumptions  of  Proposition 
o 

4.1  and  u0  e  H  (■) ,  then  there  exists  a  unique  smooth  solution  of  (S.1)  on  a  maximal 
interval  (0,Tg)  *  R,  0  <  TQ  <  ■*. 

Our  objective  is  to  study  the  problem  of  the  development  of  singularities  (shocks)  of 
the  solution  in  finite  time  su<3i  that  a  physically  meaningful  entropy  condition  will  be 
satisfied  (see  (141),  assuming  that  a  local  smooth  solution  exists. 

Let  C  e  R  and  let  u(t,x)  be  a  smooth  solution  of  (5.1).  Hs  define  the 
characteristic  through  t  of  (5.1)  to  be  the  curve  x  -  x(t,5)  in  the  t,x  plane 
specified  by  the  initial  value  problem 

' P'  (u(t,x) )  (t  >  0) 

(C  e  r)  . 

It  should  be  noted  that  the  total  derivative  of  u(t,x)  along  the  characteristic  x(t,C) 
is 


(5.2) 


dx 
dt 

x(o,o  -  e 


u(t,x(t,£))  -  ut(t,x(t,£))  +  ux(t,x(t,£)> 

-  ufc{t,x(t,?))  +  <f'  (u(t,x(t,£)))ux(t,x(t,£)) 


def 


ut(t,x(t,£))  +  v>(u(t,x(t,£)))x 


Let  u(t,£)  »  u(t,x(t,£)).  If  u  is  a  solution  of  (5.1)  its  derivative  along  the 
characteristic  x(t,C)  satisfies  the  integrodif ferential  equation 


26- 


(5.3) 


du 

dt 


a'(t 


T)<t(u(T,x(t,C)))xdT 


(t  >  0) 


l  u(0,x(0,£))  »  Up (x(0,€) )  -  Up(C)  ((<!). 

The  reader  should  note  that  (5.3)  Is  no  longer  of  convolution  type,  because  of  the  term 
♦(uCr.xft.C) )  )x  •  ♦Mu(T,x(t,5)))ux(T<x(t,C))  under  the  integral. 

Let  x(t,£)  be  the  characteristic  of  (5.1)  through  £  and  define 
v(t,5)  -  x(t,€)i  note  that  v(0,£)  >1.  The  function  v  measures  the  growth  of  the 

characteristics  with  respect  to  C.  Let  *>"(•)  *  0.  According  to  Immma  2.1  of  [14]  a 
singularity  will  develop  in  ths  solution  u  of  (5.1)  in  finite  tine,  if  it  can  be  shown 
that  there  exists  a  number  f,  0  <  f  <  •  such  that  v(f,€)  <  0.  For,  in  this  case  there 
exists  a  0  <  T  <  f  and  £1  *  ?2  e  R  such  that  xtT,^)  ■  x(T,t2>  (l*.  the 
characteristics  through  and  $2  cross  at  T),  and  u(f”,x(T,C1) )  *  u(t”,x(T,C2)) . 

This  is  the  definition  of  the  development  of  a  "shock"  at  f  in  the  solution  u  of 
(5.1).  It  is  explained  in  [14]  that  this  "shock"  solution  satisfies  the  physically 
meaningful  entropy  condition. 

We  shall  therefore  set  up  a  differential  equation  for  v(t,C) .  Let 
3U 

w(t,C)  -  (t,x(t,£)),  and  note  that  w(0,€)  ■  uj(t).  Then  using  (5.2) 

2"  dt  x£(t'5’  mJt  2"  If  *'<u<t,x(t,£>)> 


*  *"(u(t,x(t,£)))w(t,£)  . 


Thus  v  satisfies  the  initial  value  problem 

(5.4)  -jjj  *  <fim  (u(t,x(t,€) )  )w,  v(0,£)  -  1  . 

The  equation  satisfied  by  w  is  found  by  differentiating  (5.3)  with  respect  to  ?, 
obtaining  the  initial  value  problem 


(5.5) 


dw 

dt 


+ 


-v(t,5)  /  a*  (t  -  T)tr(u(T,x(t,£)))u*(T,x(t,£)) 

0 

<i'(u(T,x(t,£)))uxx(T,x(t,£))]dT,  w(0)  -  u^(£)  . 
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Our  objective  la  to  uaa  the  system  of  four  nonlinear  equations  ( 5. 2) -( S. 5)  for  the 

quantities  x(t,£),  u(t,€),  v(t,C),  w(t,C)  satisfying  the  indicated  initial  conditions  to 

establish  the  development  of  a  shock  in  the  sense  described  above.  This  problem  which  is 

under  active  study  remains  to  be  solved  in  this  generality. 

We  restrict  ourselves  to  the  special  case  ♦  =  *  in  (5.1).  By  Remarks  3.5,  3.6  the 

Cauchy  problem  (5.1)  has  a  unique  smooth  solution  u  on  [0,**)  *  II  if  a  and  •f  satisfy 

the  assumptions  of  Theorem  3.1,  0  <  a(0)  <  1,  and  if  lu  I  is  sufficiently  small. 

HZ(») 

In  the  special  case  ♦  •  *'  we  can  use  the  method  of  MacCamy  [11]  and  Dafermos  and  Nohel 
[4],  introduce  the  resolvent  kernel  k  of  a*  defined  by  the  equation 
(k)  k(t)  ♦  (s**k) (t)  -  -a*(t)  (0  <  t  <  •)  , 


(5.6) 


ut  +  V(u)x  +  k(0)u  +  k**u  -  k<t)uQ(x) 


and  write  (5.1)  in  the  equivalent  form 

(o  <  t  <  *,  x  e  i) 

u(0,x)  -  Uq ( x)  (x  e  R)  . 

Note  that  since  a  satisfies  assumptions  (a),  k(0)  «  -a'(0)  >  0.  The  method  of  Lemma  4.2 
applied  to  equation  (k)  shows  that  since  0  <  a(0)  <  1,  one  has 
(5.7)  k.k'  e  L^O,-)  . 

Remark.  (5.7)  also  holds  if  a(t)  ■  a—  ♦  A(t),  a(0)  *  1 ,  a.  >  0,  and  A  satisfies 
assumptions  (a). 

To  establish  the  development  of  singularities  in  a  smooth  solution  u  of  (5.1)  with 

>1 1  1  (equivalently  of  (5.6))  we  study  the  system  of  nonlinear  equations  corresponding  to 

(5.2)-(5.5>.  In  this  case  it  is  easily  seen  that  the  quantities  x(t,C), 

9  9 

u(t,£)  »  u(t,x(t,C>),  v(t,5)  -  ^  x(t,5),  w(t,C)  “  u(t,x(t,C) )  satisfy  the  initial 


value  problem 


(5.8) 


-  <t'  (u(t,x(t,£) ) ) 

♦  k(0)u  +  /  k' (t  -  T)u(T,x(t,S))dT  -  k( t)un (x) 
dC  o  0 

•gj-  ”  *j"(u(t,x(t,5))> 

t 

+  k(0)w  +  (/  k* (t  -  t)u  (t,x(t,5))dT  -  u' (t,x(t,? ) ) )v  -  0 


[  x(0,5)  -  C,  u<0,£)  -  uQ  (C) ,  v(0,5)  -  1,  w(0,5)  -  u^(5>  . 

In  recent  joint  work  with  Halek-Hadani  [15]  we  have  established  the  following  result* 
its  proof  which  uses  (5.7),  (5.8)  and  the  general  strategy  for  the  formation  of  shocks 
outlined  above  will  appear  elsewhere.  A  similar  result  was  stated  by  MacCamy  in  his 
lecture*  see  note  in  [6]  . 

Theorem  5.1.  Let  a  (0  <  a(0)  <  1)  satisfy  aaumptlone  (a).  In  addition,  let  the 
resolvent  kernel  k  _of_  a*  satisfy 

(5.9)  k(t)  »  0,  k'(t)  <  0  (0  <  t  <  •)  . 

Let  f  9  C2  (■) ,  *(0)  -  0,  *■(•)  >  0,  *"(•)  >  0.  tot  Uj,  e  C2(«).  If  u’(0  <  o  and 
|u^(€)l,  CeR,  is  sufficiently  large,  every  (necessarily)  smooth  solution  u(t,x)  of 
the  Cauchy  problem  (5.6)  (<■»“>  (5.1)  if  ♦  £  i f)  will  develop  a  shock  in  finite  time.  If 
•f"  ( • )  >  8  >  0 ,  an  upper  bound  for  the  time  at  which  a  shock  develops  is 


* 


_ 1_ 

k(0) 


u|,(C)8 

X°9  yi’oa)6  +  k(0)  * 


The  following  considerations  provide  examples  of  kernels  a  in  (5.1)  (<«•“>  (5.6)  if 

♦  5  ^>)  for  which  Theorem  5.1  can  be  applied. 

Remark  5.2.  If  a(t)  "8a  (0  <  8  <  1),  a  simple  calculation  shows  that 

.  -at  1-8) t 

k(t)  "  a8e  ,  and  evidently  the  inequalities  (5.9)  are  satisfied  in  the  strict 

sense . 
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More  generally,  one  has  the  following  result  established  by  elementary  consideration 


from  the  resolvent  equation  (k) . 

Lemma  5.3.  Let  a  8  C2[0,m),  (-1)^a*^(t)  >  0,  j  «  0,1,2  (0  <  t  <  "),  and  assume  that 

a1 (t)a' (0)  -  a"( t)  <0  (0  <  t  <  «)  . 

2  1 

Then  k(t)  >  0,  kMt)  <  0  (0  «  t  <  «).  If  also  0  <  a(0)  <  1,  and  a  e  W  '  [0,“),  then 

k,k'  e  l’jo,-). 


m  -at 

Corollary  5.4.  Let  a(t)  «  l  B  e  ,  8.  >  0,  a  >  0»  i £ 

j-1  3  3  3 


<*.  >  I  a.8.  U  -  1,2,... ,m)  , 

j-1  3  3 

r  1 

then  k(t)  >  0,  k*(t)  <  0  <0  <  t  <  •>»  if  also  I  8  <  1  then  k,k*  6  L  (0,-). 

j-1  3 

Finally  we  remark  that  the  general  approach  used  to  prove  Theorem  5. 1  can  also  be  used 
to  show  that  if  by  contrast,  u^(C)  >  0,  and  if  the  other  assuraptione  of  Theorem  5.1  hold, 
no  singularities  develop  in  the  solution  u  in  finite  time.  Thus  Theorem  5.1,  together 
with  this  remark,  form  the  analogue  for  the  conservation  law  with  memory  (5.1)  of  Lax's 
classical  result  for  the  conservation  law  (1.2V. 


i 


Appendix 


For  the  convenience  of  the  reader  we  state  and  prove  the  following  elementary 
inequalities  which  were  used  in  the  proof  of  Theorem  3.1  and  which  are  generally  referred 
to  as  Poincare  inequalities. 


1 

Lemma  A.  1 .  Let  g,g'  e  LZ(0,1),  g  real,  and  let  g  “  /  g(x>dx.  Then 

0 


/  g2(x)dx  <  g2  +  /  [g*(x)]2dx  , 

0  0 


In  particular,  if  g  ■  0,  then 


/  g2(x)dx  <  /  [g*(x>]2dx  . 

0  0 


Proof.  Take  0<x0<x<1.  Then 


g(x) 


x 

9<*o>  “  /  g'(E)d£  , 


and  by  Cauchy-Schwartz 


( g(x)  -  g(x0))2  <  (x  -  x0)  /  <g'(5))2d5  <  /  (g'{5)>2dC  . 

*0  ° 


Thus 


/  (g(x>  -  g(Xg))2dx  -  /  g2 ( x) dx  -  2g(x0)g  +  g2(x0) 

0  0 

1 

<  /  [g* (x) ] 2dx  . 

0 


By  the  continuity  of  g  choose  xQ  such  that  g(xQ)  “  g,  and  the  first  Inequality  Is 
immediate. 
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Lemma  A. 2.  Let  g,g’,g“  e  L2(0,1),  g  real ,  g(1)  »  g(0),  g  -  /  g(x)dx  -  0.  Then 

0 


g2(x)  +  [g'(x>]2  <  /  lg"(x)]2dx  <0  <  x  <  1)  . 
0 


Proof.  Let  0  <  y  <  x  <  1j  we  have 

x 

g'(x)  -  g'(y)  -  /  g"(£)dS 

y 

Squaring  both  sides  and  using  Cauchy-Schwartz  gives 


lg'(x)]2  *  [g'(y)]2  -  2g'(x)gMy)  <  /  [g“(x)]2dx  . 

0 


Integrating  with  respect  to  y  over  [0,1]  and  using  /  g'(y)dy  -  g(1)  -  g(0)  »  0 

0 

have 


1  1 

[g'(x)]2  +  /  lg* (y) )2dy  <  /  [g"(x)l2dx  . 
0  0 


Since  g  -  0  the  conclusion  follows  from  the  inequality 


lg< x) 1 2  <  J1  [g'(y)]2dy  , 
0 


the  proof  of  which  is  contained  in  that  of  Lemma  1. 

Application.  If  we  X(M,T)  defined  in  the  proof  of  Proposition  4.1,  then 


1 


2(t,x)  +  w2(t,x)  +  w2(t,x)  <  /  lw2x(t,x)  +  w2x(t,x)]dx  (0  <  x  <  1) 


we 
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